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Abstract 

We propose an non-standard method to calculate non-equilibrium physical 
observables. Considering the generic example of an anharmonic quantum 
oscillator, we take advantage of the fact that the commutator algebra of second 
order polynomials in creation/annihilation operators closes. We solve the 
von Neumann equation for the density-operator exactly in the mean field 
approximation and study the time evolution of the particle number and the 
expectation value {X"^). 
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INTRODUCTION 



The consistent calculation of equilibrium finite-temperature physical observables requires 
a resummation of the leading temperature contributions along the lines of the Braaten- 
Pisarski resummation program established already some years ago. The application of 
those ideas to scalar theories, QED, QCD and the electroweak model has been extremely 
successful so far, although the original resummation may have to be modified slightly near 
the light-cone for models involving massless particles. On the other hand, many interest- 
ing realistic physical scenarios like, e.g., phase transitions and particle production in a 
plasma are clearly out-of-equilibrium or only nearly-equilibrium systems in which case the 
field theoretical framework lacks a systematic basis and the straight-forward application of 
basically linear response methods gives rise to potential inconsistencies, such as pinch singu- 
larities 0. In this note we consider the toy- model of a quantum-mechanically anharmonic 
oscillator which can be solved exactly. 



THE THERMAL ANHARMONIC OSCILLATOR 

The model under consideration is characterized by the Hamiltonian 

H = H, + Hj, H,= '-P^ + ^X\ Hj = ^X\ (1) 

the quartic interaction-therm mimicking a self-interaction term as in the field-theoretical 
analog. Since there quadratic Hamiltonians are the only ones which can be integrated out 
exactly, we will also here focus on the so-called mean-field approximation. The basic idea 
is to expand the Hamiltonian operator around an expectation value that at finite tem- 
perature T (and for mixed states in general) is given by the trace of the product of the 
corresponding operator with the normalized density-operator of the state. For the Hamil- 
tonian above being parity even, the natural expansion parameter is SX"^ = X"^ — (X^) . 
Dropping the (5X^)^ term, the quadratic contribution of the interaction Hamiltonian, 
jjquadr. _ y ^)^^3 (^x'^) (2X^ — (X^)), has a form analogous to Hq with time-dependent fre- 
quency. In terms of creation and annihilation operators, a = \/2ujX + i^2/ujP and a) , the 
effective free Hamiltonian ffJ^ac^'*- under study is thus composed of 

Hf""'''- = ^f{t){a'^^ + 2a)a + a^ + 1- f{t)) and = uj{a) a + ]^) (2) 

where f(t)/uj = (X^) = Z^^Tr(pX^) has to be determined self-consistently. Note that 
diagonalizing the quadratic Hamiltonian would render the vacuum definition a |0) = time 
dependent. 

The density operator for the system in thermal equilibrium 

peq = Z^q exp{-H/T), Z^q = Tr(pe<?), (3) 

which minimizes the equilibrium entropy functional 5* ~ Tr(plogp), is a particular solution 
of the von Neumann equation of motion 
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'jtP=lH,p]. (4) 

Solutions of that equation do not uniquely determine the density-operator but admit in 
general additive static solutions that are functionals of eventual constants of motion. I want 
to point out that the equilibrium entropy functional remains constant in time for every 
solution of @. For simple, e.g. zero-dimensional integrable or infinite dimensional non- 
interacting, systems it is argued that thermodynamic information loss has to be induced 
by a second averaging procedure, e.g. coarse graining in phase space [H or in temperature 
1^, or time smoothing. For more complicated systems as quantum field-theoretical mod- 
els or classical many-body systems the second averaging process is hidden in asymptotic 
approximations, and particularly in cutting the Bogoliubov-Born-Green-Kirkwood-Yvon hi- 
erarchies 0, generally resulting in a kinetic equation. We want to mention that besides 
those somehow artificial and pragmatic methods there exists, based on a couple of evident 
assumptions, a systematic approach to clearly distinguish thermodynamic (the 'relevant') 
and microdynamic (the 'irrelevant') part of the dynamics which, however, has not yet 
been applied systematically within the framework of quantum field-theory so far. We will 
not address these open questions further but pursue the modest task to find a solution of 
the von Neumann equation. 

The biquadratic form of (0) in a, suggests the solution p 

p = exp{Ba^'^) exp{A a^a + gl) exp{Ba'^), Z = Tip, (5) 



where A, B and g are time dependent functions to be determined, and the hermiticity of the 
density operator p = p^ necessitates A, g real. On account of the basic relation 

the set of operators 5 = |l, a^a, a?, a)"^^ forms a closed commutator algebra 



a, 
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a' , a 



2{l + 2a^a) \a^^ , a] = -2a^^ \a^,a^o\=2a^ (6) 



which is the essential property that allows the evaluation of the expectation values in what 
follows. Using relations such as 

gBat2^2g-Bat2 ^ ^2 _ 35(1 + 2a) o) + AB'^o)'' 

it is elementary algebra to commutate the summands of H to the r.h.s. of (^ and to compare 
with their counterparts in the time derivative. In the resulting system of ordinary first order 
differential equations {B = x + iy), 

A = 2g = AXfy (7) 
x = 2y{l + \f) + AXfxy (8) 
-y = 2x(l + A/) + 2\fix' - y') + A(l - e'^)/2 (9) 

where the dot denotes differentiation with respect to the dimensionless time scale r = ut, 
one only needs to express f(t) in terms of A, x, y. For this purpose, it is useful to study the 
averaged elements of S and to commutate a, through p. Using the cyclicity of the trace, 
this leads to a linear system of equations for the corresponding averages with the solution 
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1 ^ (4r^ + e^-e^^) 1 , 2{x + iy) 
n = -Tiipa^a) = — , -Tr(pa ) = (10) 

where = + y"^, N = {1 - e'^f - Ar'^, thus f{t) = (1 - e^^ + Ar'^ + Ax)/{2N). 

The system (0-^ can be solved exactly. An integral of motion can be obtained by forming 
XX + yy and eliminating / by means of 

r2 = i(l + e2-4 + Ce^), (11) 
with an integration constant C. A second integral arises from eliminating y 



with a second constant of motion D. Finally, rescaling the thermal average by z = — /(C + 
2)/2, the time evolution is governed by the first order differential equation 

W). Piz) = l-^ + 2Cz-Az' + ^{z'-D) (13) 

which may be solved in terms of Jacobi integrals. Note that the solution of P{z) = 
determines the static limit of the average {X'^). 

The physical content of (|13|) may best be displayed by choosing particular initial con- 
ditions. Consider, e.g., that at t = the system be prepared in the non-self- interacting 
(A = 0) equilibrium state, to wit A{0) = Aq = —uj/T, x(0) = ?/(0) = 0, corresponding to the 
constants of motion 

n _ p-^o^2 

_ _^^e^o _L o-^o^ - S± r L _L \n _L oA)^2 

In that case 



Co = -(e^°+e-^"), Do= ' ' (8 + A(l + e-^")^j ■ (14) 



/(t) = ^^YZ^ (l - ^"""^ - ^- ^^'m\m)) (15) 
where sn is the Jacobian Elliptic Function [H with 



and module m = where 



= 11^1^(2(1 - e^°) + 3A(1 + e^°)± 



4(1 - e^o)2 + 12A(1 - e2^o) + A2(l + eA>)2 . (17) 



The expectation value f{t) performs an anharmonic oscillation with period V = K{'m)/Q 
where K{m) is the Complete Elliptic Integral of the first kind §]. We observe that also the 
occupation number corresponding to Hq varies in time. 
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no 



no + - 



l_e^o)4 
1 \"i 



2(1 - e^^°) - M_sn^(2fit|m) sn^(2fit|m) 



(18) 



which is not surprising since the particle- number operator Hq/uj no more commutes with 
the quadratic part of H. However, in the thermal average of f/''?""'^^-^ 



h 



(19) 



time dependence cancels out in a non-trivial way, as it should be since we have solved 
the von Neuman equation exactly. In fact, the two integrals of motion (^,0) are already 
sufficient to determine the energy h of the anharmonic oscillator 



C 



+ 



2XD 



2(C + 2) (C + 2)2 

for arbitrary initial conditions. 

It is interesting to study the high temperature limit where the period approaches 
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r(5) 
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2XT 



^4' 



C(T- 



(21) 



which is a valid approximation for 2AT ^ un. The period is a non-analytic function in 
the coupling constant in the high temperature limit. Physically this can be understood as 
follows. For sufficiently high T the energy levels with En ~ T are all occupied. In that 
regime the shape of the potential is dominated by the quartic part, the quadratic part being 
rather a perturbation of it. 



SUMMARY AND CONCLUSION 

We studied the dynamics of a thermal non-equilibrium anharmonic oscillator. Expanding 
the Hamiltonian around the temperature and time dependent expectation value {X'^), the 
von Neuman equation can be solved exactly for the density operator. We find that the 
system performs an anharmonic oscillation. Specifying the initial conditions to be the non- 
interacting thermal equilibrium configuration, the time evolution is studied in some detail. 

The basic ingredient in the calculation is that the commutator algebra of second order 
polynomials in creation/annihilation operators closes. This allows to commutate operators 
though the density operator which, together with the cyclicity of the trace, leads to a set of 
linear equations. It would be interesting to extend this method to go beyond the mean-field 
approximation which amounts to consider fourth order polynomials. 
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